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ABSTRACT 

We studied near-infrared disk fractions of six young clusters in the low-metallicity environments 
with [0/H] ~ —0.7 using deep JHK images with Subaru 8.2 m telescope. We found that disk 
fraction of the low-metallicity clusters declines rapidly in <1 Myr, which is much faster than the ~5- 
7Myr observed for the solar-mctallicity clusters, suggesting that disk lifetime shortens with decreasing 
mctallicity possibly with an ^10 z dependence. Since the shorter disk lifetime reduces the time 
available for planet formation, this could be one of the major reasons for the strong planet-metallicity 
correlation. Although more quantitative observational and theoretical assessments are necessary, our 
results present the first direct observational evidence that can contribute to explaining the planet- 
metallicity correlation. 

Subject headings: Galaxy: abundances - infrared: stars - open clusters and associations: general 
- planetary systems: protoplanetary disks - stars: formation - stars: pre-main 
sequence 



1. INTRODUCTION 

Because planets are formed in protoplanetary disks, 
the lifetime of protoplanetary disks (rdi s k) is thought 
to be one of the most fundamental parameters directly 
connected to the planet formation probability (P p i). In 
particular, in the current standard cor e-accretion sce- 
nario of gas giant planet formation (jSafronovl 119691 : 
ILissauer k, Stev enson 20"07l). a long disk lifetime of more 
than 10 Myr was thought to be necessary to produce the 
planetesimals that eventually be come the rocky core of 
the g iant planets ( e.g., 1 00 Myr; ILissauer! |1993|) . How- 
ever, IStrom et al.l (|1989! ) studied the time variation of 
the frequency of disk-harboring stars in Taurus molecular 
cloud that have .ftf-band excess and suggested that the 
disk lifetime is in the range from <C3 Myr to ^10 Myr, 
which is much shorter than was thought to be neces- 
sary. After their pioneering work, a notion called the 
"disk fraction", which is defined as the frequency of 
near-infrared (NIR) to mid-infrared (MIR) excess stars 
(disk-harboring stars) within a young cluster, was devel- 
oped and widely used to study the disk lifetime (e.g ., 
IIada|[T999T lHaisch et all [200lt iHernandez etHI [20071 ) . 
In all studies, the disk fraction is found to decrease as 
a function of cluster age, showing that the disk lifetime 
is about 5-7 Myr for nearby young clusters. Because the 
observation with NIR/MIR is sensitive to only heated 
dust located at the inner disk with a stellocentric dis- 
tance of ~0.1 to a fcwAU, there still remained a possi- 
bility that the disk lifetime derived with NIR/MIR ex- 
cess emission is only for the limited region of the disk. 
However, from the (sub-)millimeter continuum observa- 
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tio n of cold dust (~10K) in the outer disk (>50AU) 
by I Andrews fc Williams! (|2005[ ). the lifetime of the outer 
disk is also found to be similar, thus the entire disk is 
now thought to disappear nearly simultaneously. After 
the discovery that Tdisk is less than 10 Myr, it is widely 
accepted that giant planet cores must form in this short 
time for any theoretical models (|Lissauer &: Stevenson! 
120071) . 

Altho ugh as many as over 450 exoplanets are now 
known (jSchneider et al.1 l2009T l . these planets are unex- 
pectedly found to have a lot of diversity in terms of 
e.g., mass, orbital period, and eccentricity. Perhaps the 
most telling discovery in exoplanet study is the "planet- 
metallicity correlation", the higher probability of a star 
hosting a giant pla net with increasing metallicity (e.g., 
iFischer et al.1 120051 ). suggesting that mctallicity could 
be the most crucial parameter for giant planet forma- 
tion. Although this correlation is generally interpreted 
as a natural consequence of the core-accretion scenario, 
the origin of it is still not well-explained despite de- 
tailed analyses by either theoretical studies assuming so- 
lar m etallicity or observations of solar met allicity regions 
(e.g. JIda fc Lm1l2004bl : IWvatt et al.ll2007T) . Therefore, it 
is necessary to directly study protoplanetary disks under 
different mctallicity conditions to find any clues of the 
planet-metallicity correlation and to test theories of it. 

Since our Galaxy is known to have a metallicity gra- 
dient with lower mctallicity for larger Galactic radius, 
rangi ng from ~— 1 dex to ~+0.5dex (e.g.. !Rudolph et all 
f2006h . we can explore different metallicity regions at 
large distances with high-sensitivity observations with 
large telescopes. As a first step, we focused on the low- 
metallicity environment in the outer region of our Galaxy 
{R g > 15 kpc) to study a region with significantly lower 
mctallicity ([O/H] ~ — ldex). As a first result, we de- 
rived the disk fra ction of two clus ters in Digcl Cloud 
2 at R g = 19 kpc (|Yasui et alj 120091 hereafter Paper I), 
which are low metallicity ([O/H] ~ —0.7) and very young 
(~0.5 Myr old) clusters. Disk fractions for both clusters 
were found to be quite low in spite of their very young 
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age, suggesting that the disk lifetime in low metallicity 
is short. In this paper, we summarize the disk lifetime in 
low-metallicity environments with our all targets in the 
outer Galaxy and discuss the implication for the planct- 
metallicity correlation. 

2. TARGET SELECTION AND OBSERVATION 

We searched the literature for star-forming regions in 
the outer Galaxy and with [O/H] < —0.5. Among the 
listed 10 candidate star-forming regions, we selected four 
regions at R g > 15kpc and in the second quadrant of the 
Galaxy as an initial set of samples, mostly by consider- 
ing the visibility from the Northern hemisphere. We also 
checked whether the candidate regions have at least one 
associated stellar aggregate using Two Micron All Sky 
Survey (2MASS) data. As a result, we selected six clus- 
ters (Table ED: on e cluster in each Sh 2-207 a nd Sh 2-208 
dBica et al.fl2003f ) . and two in each Sh 2-209 (IKlein et al.l 
[2005D and Digel Cloud 2 (IKobavashi etaD 12008ft . The 
average metallicity of all target clusters is ~— 0.7dex, 
which is significantly lower than the solar metallicity. 

Deep JHKs-b&nd images of the target star-forming 
regions were obtained with the 8.2 m Subaru telescope 
equipped with MO IRCS (Multi-Objec t InfraRed Camera 
and Spectrograph; iSuzuki et al.ll2008h , which has a wide 
field of vi ew (4' x 7') and uses t he MKO NIR photomet- 
ric filters (|Tokunaga et a l. 2002). The weather condition 
was photometric with excellent seeing (~0'.'3-0'.'4) for the 
Cloud 2 and Sh 2-209 observations on 2006 September 
3 UT, while a little cirrus was present with seeing of 
-1" for the Sh 2-207 and Sh 2-208 observations on 2006 
November 8, 2007 November 23, and 2008 January 14 
UT. Total integration time for each wavelength band was 
~500-1000sec. After the raw data were reduced with 
standard IRAF procedures, JHK aperture photometry 
was performed. For most clusters aperture photometry 
with IRAF apphot was enough even for such large dis- 
tances because stellar images arc sufficiently smaller than 
the stellar separations. Only for two dense clusters in 
Sh 2-209, we performed photometry with point-spread 
function fitting using IRAF daophot with scripts for au- 
tomated photometry, "autodao"H 

Photometric standard stars that were obtained at the 
similar airmass as the object fields are used except for 
Sh 2-207 and Sh 2-208 clusters, for which 2MASS stars 
in the fields are used after converting the 2MASS mag- 
nitudes to the MKO mag n itude s with the color transfor- 
mations in iLeggett et al.l ()2006l ) . The achieved limiting 
magnitudes (lOer), J ~ 20-22 mag, H ~ 19-21 mag, and 
Ks — 18.5-21 mag, correspond to <0.1M S . This mass 
detection limit, sufficiently less than 1 M , enables a di- 
rect comparison of disk fractions in such distant clusters 
with tho se of star-forming c lusters in the solar neighbor- 
hood (cf. lHaisch et al.ll2001l ). For each embedded cluster, 
we identified clus ter mem bers in the same way as for the 
Cloud 2 clusters (jPaper J ): the extinction criterion {Ay) 
for identifying cluster members was determined for each 
cluster considering the Ay distributions in the cluster 
region and the field region (see details in lPaperll ). The 
details for each cluster will be published in separated pa- 
pers. 

6 Copyright (C) 2008-2009 Noriyuki Matsunaga 
(http: //www. ioa. s .u-tokyo . ac . jp/~nmatsuna/ Japanese/sof tware/aut 
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Fig. 1. — Disk fraction as a function of cluster age. JHK disk 
fractions of the young clusters with low metallicity arc shown by red 
filled circles, while those of young clusters with solar metallicity are 
shown by black filled circles. The typical systematic error of ages 
for the solar metallicity cluster is shown at the top-right corner (see 
the main text for explanation of the gray horizontal lines for three 
of the clusters). The black line shows the disk fraction evolution 
under solar metallicity, while the red arrow shows the proposed 
JHK disk fraction evolution in low-metallicity environments. Note 
that these lines are drawn by eye and are not derived with rigorous 
fitting. 



3. DERIVATION OF AGE AND DISK FRACTION 

3.1. Clusters in Low Metallicity 

Following the method bv lMuench et al.l (|2000l ). we first 
derived the ages of clusters by comparing observed the 
-ftT-band luminosity function (KLF) with model KLFs of 
various ages constructed from the canonical initial mass 
function (IMF), mass-luminosity (M-L) relation and dis- 
tance of the cluster. The detail of the model fitting is 
described in lPaper 1 in the case of Cloud 2. We applied 
the same procedure to all clusters. We confirmed that 
the variation of the IMF and M-L relation in the rele- 
vant metallicity range is small enough so that the ages 
of the clust ers can be estima ted within the uncertainty 
of <lMyr (| Yasui et all 120061 ). which is sufficient for the 
present study. Estimated ages of most targets are found 
to be very young, ^lMyr, but one cluster, Sh 2-207, 
shows a relatively old age (2-3 Myr), as expected from 
the 2MASS image, which shows a less-extincted cluster. 

Next, we used JHK color-color diagrams to derive 
the disk fractions for each cluster. The method we used 
to determine cluster infrared excess fr actions f rom color- 
color diagrams is described in detail in lPaperll along with 
extensive discussions on the uncertainties. In Table [1] 
we list the derived disk fractions and ages for the entire 
sample clusters, and plotted the results in Figure [T] The 
uncertainties of the disk fraction reflect Poisson errors. 
Only for the Sh 2-207 cluster we derived the disk frac- 
tion by considering the contamination from the field stars 
because i t is difficult to identify cluster members with 
extinction of this older cluster. The number 
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TABLE 1 

JHK DISK FRACTIONS OF EMBEDDED CLUSTERS IN LOW-METALLICITY ENVIRONMENTS. 



Cluster 


log [0/H] a 


Rg /D b 


Age 


Disk Fraction 


M lim c 




(dex) 


(kpc) 


(Myr) 


(%) 


(M e ) 


Cloud 2-N 


-0.7 (1) 


19/12 (i) 


0.5-1 


9±4 (5/52) 


0.06 


Cloud 2-S 


-0.7(1) 


19/12 (i) 


0.5-1 


27±7 (16/59) 


0.06 


Sh 2-207 


-0.7 (2) 


12/4 d (ii) 


2-3 


5.1±4.6 (1.2/23.4) e 


0.07-0.08 


Sh 2-208 


-0.8 (2) 


12/4 d (ii) 


0.5 


19±5.5 (12/63) 


0.05 


Sh 2-209 main-f 


-0.6 (2) 


17.5/10 (iii) 


0.5-1 


10±0.8 (163/1605) 


0.08-0.09 


Sh 2-209 sub-f 


-0.6 (2) 


17.5/10 (iii) 


0.5-1 


7.1±1.2 (35/494) 


0.08-0.09 



References. — CD fLubowich et ail l f20"oll ), (2") |Caplan et al.| J2000f) ; fn iKobayashi et al.| J2008I) , qp lFich et"aT1 il990l) . (ni) IChini Wink! fa 984) . 



Note. — "Metallicity of the observed regions. The solar metallicity is assumed as log(0/H)Q — 8.7 dAsplund et al1|2009D . b Adopted distance 
(D) and Galactic distance (R g ). We assumed the Galactic distance of the sun as Rq = 8.0 kpc (lReidlll993h . c Mass detection limit of the data. 

For both Sh 2-207 and Sh 2-208 clusters, wc adopted the kinematic distances because model KLFs with the photometric distances (D 8— 9 kpc, 
R g ~ 15.5— 16.5 kpc) do not match the observed KLF at all, while those with the kinematic distances (D 4 kpc) match very well. e See detail in 
the main text. -^Sh 2-209 has two clusters, which are located near the two peaks of millimeter continuum (Klein ct al. 2005): hereafter wc call the 
larger cluster at the northeast side and the smaller cluster at the southwest side as "main cluster" and "sub-cluster" , respectively. 



TABLE 2 

JHK DISK FRACTIONS OF EMBEDDED CLUSTERS WITH AGES OF >5 MYR IN 
THE SOLAR. NEIGHBORHOOD. 



Cluster 


Agc a 


Disk Fraction 


M lim " 


Filter c 


Rcf d 




(Myr) 


(%) 


CM©) 






Orion OBlb 


5 (1) 


8.9±3 


0.1 


2MASS 


1 


Upper Sco 


5(2) 


32±4 


0.1 


2MASS 


3 


•q Cham 


6(4) 


28±12 


0.08 


2MASS 


5 


Orion OBla 


8.5 (1) 


3.6±2.5 


0.1 


2MASS 


1 


NGC 7160 


10 (5) 


6.3±3.6 


0.4 


2MASS 


6 



Note. — "References for the ages are shown in the parenthesis. 
Although there are two age estimates for Orion OBlb and OBla by 
using different isochrone models, wc show the average value here. 
'Mass detection limit of the data. c The photometric system of the 
obtained JHK data. ^References for the JHK photometric data. 
For Upper Scorpius and NGC 7160, we used JHK magnitudes in 
2MASS All-Sky Point Source Catalog for the stars listed in the 
references. 

Ref ere nces. 
(2002), (3) Carpenter' ct al. (2006), (4) 



(1) IBriceno et al.l <2005l). (2) IPreibisch et all 
r ct al. (2006), (4) Luliman & Stceghs (2004), 



(5) Sicilia-Aguilar ct al. (2005), (6) S icilia-Ag uilar ct aE] J2006I ) 



of stars in the Sh 2-207 cluster region is 38 and of these 
only 2 have disks. The normalized number of stars in the 
field region is estimated to be 14.6±0.9 with that lying in 
the infrared excess region of the diagrams to be 0.8±0.2. 
Therefore, the final disk fraction of the Sh 2-207 cluster 
is estimated at 5.1% ± 4.6% ((2 - 0.8)/(38 - 14.6)). 



7160; see Table EJ. 

Although the ages of these clusters are derived with 
other methods than KLF fitting (e.g., with spectroscopy 
and imaging data using the H-R diagram), we derived 
the ages of the clusters with KLF fitting for better con- 
sistency. Because the KLF fitting is reliable only for 
good signal-to-noise imaging data which cover the peak 
magnitude of KLF, we derived ages of Trapezium, Up- 
per Sco, and NGC 7160 clusters, since these are the only 
clusters for which we could find such data in the liter- 
ature. The resultant ages are identical to those in the 
literature in spite of larger uncertainties with increasing 
ages as shown with thick gray horizontal lines in Fig- 
ure 1. Considering that ages of reference clusters esti- 
mated with various pre-main-se quence models have sys- 
tematic uncertainty of ^lMyr l)Haisch et al.ll200ll ). the 
age estimate of clusters with low-metallicity with KLF 
fitting in § 13. II appears to be reliable. 

For all clusters, the JHK disk fractions were estimated 
in the same way as for the clusters in low-metallicity en- 
vironments. The resultant disk fractions and ages for 
solar-metallicity clusters are shown in Figure 1 with black 
circles. They show decreasing disk fraction with increas- 
ing age: the decline is rapid up to 5-7 Myr and stays 
almost flat beyond this point with disk fraction slowly 
declining from 10% to 5%. These characteristics are to- 
tally identical to the latest MIR disk fractions based on 
ground-based L-band observations and Spitz er observa- 
tions (see Figure 14 in lHernandez et aT1l2007[) . 



3.2. Reference Clusters in Solar Metallicity 

For comparison, wc also derived the JHK disk frac- 
tion and age of various embedded clusters in the solar 
neighborhood in a uniform way using the photometry 
data in the literature. We chose publications with the 
following criteria: (1) JHK photometry data set of clus- 
ter members in the same filter system are available for 
a reliab le estimate and (2 ) the mass detection limit is 
<1 M (jHaisch et al.ll200lD to sufficiently cover the most 
typical mass range in order to have the most representa- 
tive disk fraction for each cluster. In addition to seven 
yo ung (< 5 Myr) clusters in the solar neighborhood used 
in IPaper 11 (see Table [5]), we derived disk fractions for 
additional five clusters with older ages (>5Myr; Orion 
OBlb, Upper Scorpius, n Cham, Orion OBla, and NGC 



4. DISK LIFETIME IN LOW-METALLICITY 
ENVIRONMENTS 

Figure 1 clearly shows that the disk fractions of all 
the clusters in low-metallicity environments are less than 
30%, which is significantly lower than nearby embedded 
clusters of similar age. Also, the disk fraction with low 
metallicity seems to reach the first lowest level (~10%) 
at the age of ~1 Myr, which corresponds to the disk frac- 
tion with solar metallicity at ~5-7Myr. In combination 
with the fact that there are no embedded clusters of so- 
lar metallicity whose disk fractions are this low at the 
same ages, our results show that disk fraction strongly 
de pends on metallicity, confirming our earlier suggestion 
in lPaper"H . 

Following the detailed discussion in IPaper J . our re- 
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suits suggest that the entire disk dispersal occurs in the 
very early phase in low-metallicity environments, while 
initially (at age = OMyr) an inner disk as well as outer 
disk exists even under low mctallicity. The disk lifetime 
in low metallicity is estimated at ~1 Myr from Figure 1, 
while that in solar metallicity is ~5-7Myr, suggesting 
that very short disk lifetime in low-metallicity environ- 
ment. Assuming that the disk lifetime-metallicity rela- 
tion is described as Tdisk oc 10 aZ , the constant a is esti- 
mated to be 1 ± 0.5 despite that we have only two mea- 
surement points (Odcx and — 0.7dex) at present. The 
uncertainty of the constant comes from the uncertainties 
of metallicity and age of each cluster. Note that our NIR 
observations are sensitive to only heated dust in the in- 
ner disk (see detail in § 1) and millimeter /sub- millimeter 
observations of the outer disk (e.g., with ALMA) in the 
near future are very important to conclude that the above 
estimated disk lifetime is actually for the entire disk. 

The metallicity dependence of the disk lifetime may 
impose constraints on the mechanism of disk disper- 
sal, which is still under debate even for solar metal- 
licity environments. Among five major mechanisms 
of disk dispersal: planet formation, stellar encounter, 
strippi ng by winds, mass acc retion, and photoevapo- 
ration (jHollenbach et all 120001 ) . only photoevaporation 
by the radiation of the central star appears to cause 
the strong metallicity dependen ce (see discussion in 
iPaper J : lErcolano fc Clarke! l2010T l if the dominant ra- 
diati on is in far-UV (FUV, 6 eV < hv < 13.6 eV; 
e.g., IGorti fc Ho lenbachl [2009ah and /or X -rays as sug- 
geste d in IPaper J (hv > 0.1 keV; e.g.. lErcolano fc Clarke! 
I201C0 . However, because the theory of photoevapora- 
tion is still under development and with numerous pa- 
rameters, further studies are needed. The metallicity 
dependencies of disk lifetime from theories are ocl0 3Z 
and ( xlQ - 5 - 2 for FUV photoevap oration from the theo- 
ries of lGorti fc Hollenbachl (12009a!) and X-ray ph otoevap- 
oration from those of lEr colano fc Clarke! ([2010h . respec- 
tively. This is not enough to explain our results, although 
it is marginally consistent within the uncertainty (10 aZ , 
a = 1 ± 0.5). Note that s omething specifi c to the outer 
Galaxy environment (e.g.. iHavwoodl 120081 ) might be the 
reason for the rapid disk dispersal. However, considering 
that Sh 2-207 and Sh 2-208 clusters, which are not found 



to be located in the outer Galaxy region (R g = 12kpc, 
see note in Table 1), also show low disk fraction and thus 
short disk lifetime, mctallicity could be the dominant fac- 
tor. 



5. IMPLICATION TO THE PLANET-METALLICITY 
CORRELATION 

The observed planet-metallicity correlation is known 
to have strong m etallicity dependence (P p i cx 10 2Z ; 
iFischer et al.l [20051 ). Although the core-accretion model 
qualitativel y explains t he depe ndence, e.g., P p i oc 10 z in 
the case of llda fc Linl (|2004bl ). which is a deterministic 
model based on the core accretion, it is still not enough to 
expla in the steep 10 2Z dependence (see e.g. JWvatt et al.l 
2007). The short disk lifetime under low-metallicity 
suggests that planet formation becomes more difficult 
with decreasing metallicity, and therefore it may explain 
a part o f the strong planet-metallicity correlation. In 
IPaper fl . we proposed that by adding the contribution 
from the disk lifetime, most of the observed metallic- 
ity dependence of the planet formatio n probability could 
be rea sonably explained. Recently, lErcolano fc Clarke! 
(|2010f ) also estimated the correlation between P p \ and 
''disk by combining their ph otoevapora t ion mo del and 
the core-accretion model by llda fc Lml (|2004bf ) to sug- 
gest that the metallicity dependence of photoevaporation 
only plays a secondary role with an effect of < 10 02Z 
on the planet formation probability. However, there are 
more ambiguities and complexities in the planet forma- 
tion process than the the oretical planet for mation sce- 
narios now available (e.g.. llda fc Linl r2004a| ) . Also, the 
derived metallicity dependence of the disk lifetime still 
has a large uncertainty (see § [4|, and needs to be con- 
firmed by observations of clusters with higher than solar 
mctallicity. The mctallicity dependence of the disk life- 
time suggested here is the only direct observational ev- 
idence for the mctallicity dependence of the planet for- 
mation and could be an important clue to understanding 
the planet formation mechanism. 



C.Y. has been supported by the Japan Society for the 
Promotion of Science (JSPS). 
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